DIGITAL DOUBLE BARRIER OPTIONS: SEVERAL BARRIER 
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Abstract. We determine the price of digital double barrier options with an 
arbitrary number of barrier periods in the Black-Scholes model. This means 
that the barriers are active during some time intervals, but are switched off 
in between. As an application, we calculate the value of a structure floor for 
structured notes whose individual coupons are digital double barrier options. 
This value can also be approximated by the price of a corridor put. 



1. Introduction 

We consider digital double barrier options with an arbitrary number of barrier 
periods. This means that the holder receives the payoff only if the underlying stays 
between the two barriers in certain specified time intervals. While such contracts 
might make sense by themselves (as a weather or energy derivative with seasonal 
barriers, say), our motivation is to use them for the pricing of certain structured 
notes with several coupons. Such trades often feature an aggregate floor at the final 
coupon date, which increases the total payoff to a guaranteed amount if the sum of 
the coupons is less than this amount. Pricing this terminal premium requires the 
law of the sum of the coupons, which can be recovered from its moments. If the 
individual coupons of the note are digital barrier options, then these moments can 
be computed from the prices of options of the kind described above, where the sets 
of barrier periods are subsets of the coupon periods of the note. 

Recall that Monte Carlo pricing of barrier contracts is tricky, because the dis- 
cretization produces a downward bias for the barrier hitting probability. For single 
barrier options, this difficulty can be overcome using the explicit law of the maxi- 
mum of the Brownian bridge JTJ [3J. For double barrier options, the exit probability 
of the Brownian bridge is not known; see Baldi et al. [2] for an approximate ap- 
proach using sample path large deviations. These numerical challenges led us to 
investigate exact valuation formulas. 

The paper is structured as follows. In Section [2] we define the payoffs we are 
interested in and price them for a single barrier period. Section [3] extends the 
result to arbitrarily many periods of active barriers. Our main application, namely 
the pricing of structure floors, is presented in Section [4] Since our exact pricing 
formula is fairly involved, we consider an asymptotic approximation for a large 
number of periods in Section [5] 
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2. Preliminaries and pricing for one period 
We assume that the underlying (St)t>o nas the risk-neutral dynamics 

dS t /St = rdt + adW t 

with constant interest rate r > 0, volatility a > and a standard Brownian mo- 
tion W. Consider a digital barrier option with two barriers B\ OVJ and B up that are 
activated at time T > and stay active for a time period of length P > 0. At 
maturity To + P, the payoff is one unit of currency if the underlying has stayed 
between the two barriers: 



(!) ^1 : - 1 {-B low <S t <B up , te[T ,T +P]}- 

Let us denote the price of this "one-period double barrier digital" by 

(2) BD(S t ,t;{T },P,B low ,B up ,r) := e^^E^i], 

where E is the expectation w.r.t. the pricing measure P. In the terminology of 
Hui [7J, this is a rear-end barrier option, because the two barriers are alive only 
towards the end of the contract, namely between To and maturity To + P. Hui [7 
has determined the price for a barrier call of this kind. The digital case is a simple 
modification, but we go through it to prepare the calculation of the price for several 
barrier periods (see Section |3j. The value function 

f(S, t) := BD(S, t; {T }, P, B lovn B up , r) 

satisfies the Black-Scholes PDE 

9/1 2 Q 2d 2 f M q df 

with the terminal condition f(S, Tq+P) = 1, for S G (-Blow, ^up), and the boundary 
conditions f(B\ ovl ,t) = f(B up ,t) = for t e [T ,T + P]. We use the standard 
transformation f(S,t) = e ax+ ^ T U(x,T), where 

(3) x:=\og(S/B low ), t:= ±<j 2 (T Q + P - t) , 

1/2 \ n 2r 2 

o:=--^r-lj, P^-^-a , 

to transform the Black-Scholes PDE into the heat equation 

d 2 U dU 

(4) 



dx 2 dr 



The time points (0,To,T) + P) are thus converted to (|<r 2 (To + P),p, 0), where 
p := \cf 2 P is the barrier period length in the new time scale. The boundary 
conditions in the new coordinates are 



(5) U(0,t) = U(L,t)=0, re[0,p], 

where L := log(-B up /.Bi ow ). The terminal condition translates to the initial condi- 
tion 



(6) 



U{x,0) = e- ax , ie(0,L). 
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Proposition 1. For < t < To, the price of a barrier digital with barrier period 
[To, To + P] and payoff C\ at T> + P (see (fTj) ^) is 



(7) 



BD(S, t; { T 0} , P, B low , B up , r) = V^(—j g k ' ^ - ^ < 



s Ysr^i i-(-i) fc e 



/c„-qL 



t/2(t- P ) . ( kn 



(x + yy/^r-p))^ e-v 2 ' 2 dy. 



Proof. We have to solve the problem Q-([6]). First consider the rectangle (0, L) x 
(0,p). There the solution can be found by separation of variables [5J Section 4.1]: 

QQ / 7, \ 

(8) %r) = ^ fcS h e"(T)\ (x, r) G (0, L) x (0,p), 

fe=i ^ ' 

where 

' e 1 sm — ari dx\ = 2kn- 



are the Fourier coefficients of the boundary function U(x,0) = e~ ax . At r = p, 
the solution is given by Q for < x < L and vanishes otherwise. Inserting r = p 
into p| yields 



(9) U{x,p) 



= feiS^ ^K sin(^)e-(^, 0<.T<i 
1 0, x < or x > L. 



Now we solve for U in the region K x (p, ^a 2 (To + P)). There are no boundary 
conditions here, since the barriers are not active in the interval (0,To) (in the 
original time scale). The solution is found by convolving the initial condition Q 
with the heat kernel [5J 2.3.1.b]: 

1 f°° 

U(x,t) = -= U(x + yy/2(T-p),p)e-y 2 / 2 dy 
v 2ir J -co 



(10) = -j= I V ' iT - p> U(x + y^2(^p),p)e-y l 2 dy. 



Inserting Q and rearranging yields □ 

3. Double barrier digitals with arbitrarily many periods 
For n tenor dates 

< T < • • • < T„_! 

and a fixed period length P > 0, we consider a contract that pays one unit of 
currency at time T„_i + P, if the underlying has remained between the two bar- 
riers Slow and B up during each of the time intervals [Tj,Tj + P], i = 0, . . . , n — 1. 
By the risk-neutral pricing formula, the price of this "multi-period double barrier 
digital" is given by 



(11) BD{S u t;{T Q ,...,T n _ 1 },P,B low ,B^,r) := e -'( r »-*+^E 
where 



i=l 



Ci := 1 



{B Iow <S t <B up , te[T 1 _i,T 4 _i+,p]}- 
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To calculate the price, we use the coordinate change ^ again. The n barrier 
periods [Ti,Ti + P] arc mapped to [t^t,; +p], where 

n := \<j 2 (T n -i - Tj_i), i = n, 

are the images of the barrier period endpoints under the coordinate change (see 
Figure nl). Define the following auxiliary functions: 



L 



u = o 




U = 




Fourier se- 


convolution 


Fourier series, 


convolution 


ries, coef- 


Of U(;T n + 


coefficients 


of value at 


ficients by 


p) with heat 


by integrating 


T„-l + p 


integrating 


kernel 


value at r n _i 


with heat 


U(;0) 






kernel 


P T n 

u = o 




+ P 



Figure 1. Solving the boundary value problem for an arbitrary 
number of barrier periods. 



hj(ki,. . . , kj+nxt, . . . , Xj+i; yi, . . . , yj+i',x, t) 
:= ^=e"^+ l/2 1 



and 



(Vi- 



^/2(r-(T„_ J -+p))'^/2(r-(r„_ i +p))_ 

i;xi, . . . ,x J+1 ;yi,. . . ,yj;x + y j+1 ^2(r - (r„_j +p)),r„_ J + p) 



(fci , . . . , fcj+i ; xx, . . . , x j+1 ;y 1 ,...,y j ;x,T) 

. 3+ sin — e 

hj—i [h\ , . . . , fcj ; x\ , . . . , xj ; t/i , . . . , yj] , T n —j ) , 



:= | sin s i n k i+i™ c -(fc 7+ i7r/L) 2 (r-T n _ 1 ) 



with the recursion starting at 

(12) ffo(fei;»i;;a;,T") 



= le""* 1 sin sin fcii™ e -(fciT/-E")\ 



The following theorem contains our pricing formula. The first formula ( 13 ) is for 
time points inside a barrier period, whereas the second expression ( 14 ) holds for 
valuation times where the barriers are not active. 



Theorem 2. The value function (11) equals e ax+ ^ T U{x 1 r), where for < j < n, 
T n~j < t < T n~j + p, < x < L, we have 



(13) U(x,t) 



-oo «/o ^0 



-v — 



£••• E 

A.'i— fei+i =0 



3 i+i 

(fci , . . . , fcj+i ; xi, . . . , x j+1 ; yi, . . . , yf, x, r)dxi . . . dxj+xdyt . . . dy , 
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whereas for < j < n, T n -j + p < t < T n _(j_|_i) fwit/i tq :— oo), x € we fcaue 



(14) £/(z,t) 



L oo 



E- E 



i+i J'+i 

hj(ki, . . . , fcj+i;a;i,. . . ,x j+ i;yi, . . . ,y 1+ i;x,T)dxi . ..dx j+ idyi . ..dy. J+1 . 

Proof. The idea is to iterate the argument of Proposition [I] (see Figure [lj . We use 
separation of variables in the barrier periods, and convolution with the heat kernel 
for the periods in between. The required initial condition at the left boundary 
comes from the previous step of the iteration (for j — also from the payoff, of 
course) . 

For j = 0, formula ([13]) is identical to To show (fli] ) for j = 0, let p < r < 
r n _i (recall that r„ — 0) and ieR, and use (10 1 and (181) to obtain 



U(x,t) 



1 



/27T J- 
1 



%/2(t-p) 
O 

lr 



U{x + Vl V2(r-p), p)e- y "' 2 d yi 



x L—x 
' y/2(T-p) ' y/2(T-p) 



t 2 

/ E 9o{h;xi; ;x + y 1 y/2(r - p) 1 p)e~ Vl/2 dxidy 1 
J ° fc 1= o 

/oo pL °° 
/ E h o(ki;x 1 ;y 1 ;x,T)dx 1 dy 1 . 
-oo J r, n 



fcl=0 



This is (|14|) for j = 0. 

Next consider a rectangle 

(15) (r, x) E ( + p)x(0,L), 



1 < j < n. 



At the left boundary, the solution is i-> U(xj + i,T n -j). By the induction 

hypothesis, it equals (14) with j replaced by j — 1: 



(16) U(x j+ i,T n -j) = 



J_! OO oo 

E-E 

j 3 

kj\xi, ...,Xj;yi,..., yj;x j+1 ,T n -j)dxi . . . dxjdyx . . . dy 3 



The solution in the rectangle (|15|) is thus obtained by separation of variables as 

5 _ ( ^ )2(T _ Tr> _. ) 



(17) 

where 
(18) 



U(x ) r) = £ 6 fcj+1 sin(^ 



kj + i—0 

-f 

LJo 



U(xj+i,T n -j)sin 



kj+lTT 

L ' 



i I dxj- 



denote now the Fourier coefficients of x j +1 H t U(xj + i,T n -j). Inserting (16l into (18) 
and then (18) into (17) yields (13), by the definition of gj. 
Finally, consider a strip 



(19) 



(t,x) e (r n _j +p,t„_ (3+1) ) x 



1 < j < n. 
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At the left boundary, we use ( 13 ) as induction hypothesis. The solution thus 

- p and x £ (0, L) it is 



vanishes for 



(0, L), and for r = T n -j 



(20) U(x, Tn-j + p) 



E 



9j(h 



j 3 + 1 

, kj + i; xi, x j+ i;yi,. . . , y 3 \ x, r„_j + p)cfa;i 



d,Xj + idyi . ..dy 3 



As above, the solution in the strip (191 is found by convolution with the heat kernel: 
1 



U(x,t) 



%/2tt J- 



' y/2{T-{r n - j -p)) ' y/2(r-(7 



-P)) 



(v. 



3+1 J 



U(x + y j+1 ^2(r - (r n _j - j?)), r„_j - p)e v ^+ l/2 dy J+1 . 
Now insert (20), with a; replaced by a; + yj + i^2(r — (j n -j +p)), and use the defi- 



nition of hi to conclude (114 



□ 

Note that Proposition [T] corresponds to (14) for j — 0. As seen there, the 

L — 

integral J Q dxi can be done in closed form. We have not included this evaluation 
in Theorem [2] to increase its readability. 

If a different option (a call, say) with the same barrier conditions is to be priced 
instead of a digital payoff, the quantity e~ axi in (12) should be replaced by the 
appropriate payoff U(xi,0). 

4. Structure floors 

In this section we assume that our tenor structure satisfies Ti__i + P = Ti for 
1 < i < n, and define T n := T n -\ + P. We consider a structured note with n 
coupons, where the i-th coupon consists of a payment of 

( 21 ) c i = 1 {Si ow <s t <s up , te[Ti-i,Ti]}, l<i<n, 

at time Ti . These coupons can be priced by Proposition [l] (replace To by TV 
addition, the holder receives the terminal premium 

/ \ + 



In 



(22) 



F 



Ci 



at T„, where F > 0. This means that the aggregate payoff A := J^ILi of the 
note is floored at F, which is a popular feature of structured notes. While the 
individual coupons are straightforward to valuate, it is less obvious how to get a 
handle on the law of A. We now show that this law is linked to barrier options 
with several barrier periods. Indeed, the following result is based on the fact that 
the moments 



(23) 



E[A V ] = V i v P[A = i], l<v <n, 



=o 



of A are linear combinations of multi-period double barrier option prices, with 
coefficients 



(24) 



c(v, J) 



E 

supp(i) = J 



*1 } * * * t *7 



JC {l,...,n}. 



(The notation supp(i) = J means that J is the set of indices such that the corre- 
sponding components of the vector i = (ii, . . . , i n ) are non-zero.) 
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Theorem 3. The price of the structure floor (22 1 at time t < Tq can be expressed 
as 



(25) 



- r ( T "-*)E[(F-A)+] 



nA L-FJ 

-r(T„-t) £ (F _ i)P[A = i], 
i=0 



zw/iere 

(26) P[A = n] = PP(S t ,*;{To},T„-T ,Pi ow ,P U p,0). 

T/ie other point masses P[A — - i] in (|25[) can 6e recovered from the moments of A 



by solving ( |23[ ) (including v = 0, of course). The moments in turn can be computed 
from barrier digital prices by (1 < v < n) 

(27) E[A V ] = Yl c ^< J ) ' BD ( S t^ {Tj ■ ] € J}, P, Blow, P up , 0), 

JC{l,...,n} 



where the coefficients c(v,J) are defined in (24 1 



Proof. The expression (25) is clear. The event in (26) means that all of the n 



coupons (21 1 are paid. By our assumption that Tj = Tj_i + P, its risk-neutral 



probability is the (undiscounted) price of a double barrier digital with one barrier 
period [T ,T n ], which yields ([26]). To prove (27), we calculate 



V[A U ] = E 




Now observe that YijeJ^j is the payoff of a double barrier digital with barrier 
periods [Tj,Tj + P] for j e J. □ 



When calculating the value BD in (27) for, say, J = {1,2,4,5,6}, the adjacent 
barrier periods should be concatenated: Do not compute the price for five barrier 
periods of length P, but rather for two periods with lengths 2P and 3P. We did 
not include this obvious extension (barrier periods of variable length) in Theorem [2] 
in order not to complicate an already heavy notation. 

5. Approximation by a corridor option 



Theorems [2] and [3] express the price of the structure floor ( |22[ ) by iterated sums 
and integrals. Due to the factors of order e~ k i , the infinite series Ylk ma y ^ e 
truncated after just a few terms. Still, numerical quadrature may be too involved 
for a large number of coupons, so we present an approximation. Let us fix a 
maturity T = T n and assume that the n coupon periods 

T™ := [*=iT, -T], 1< % < n, 

have length T/n. For large n, the proportion of intervals during which the under- 
lying stays inside the barrier interval 

B := [Plow, P U p] 

is similar to the proportion of time that the underlying spends inside B, i.e., the 
occupation time. This is made precise in the following result, which holds not 
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only for the Black-Scholes model, but for virtually any continuous model. Note 
that the level sets of geometric Brownian motion have a.s. measure zero (cf. [H 
Theorem 2.9.6]). 

Theorem 4. Let {St)t>o be a continuous stochastic process such that for each real c 
the level set {t : St — c} has a.s. Lebesgue measure zero. Then we have a.s. 

1 - f T 
lim -Vl {StGB vter"} = / lB(St)dt. 
™^°° n i~i Jo 

Proof. For 1 < i < n, define processes (X ni (t)) a < t <T by 
X n i(t) : 



1 if t e T~ n and S u eBVueTf 
otherwise. 



Put X n := Yli=iXni- We claim that, a.s., the function X n (-) converges pointwise 
on the set [0,T]\{f : S t = B iow or S t = B up }, with limit 1 B (S.). Indeed, if t G [0,T] 
is such that St 4- then X n (t) = for all n. If, on the other hand, St G int(B), 
then t has a neighborhood V such that S u E B for all u G V, by continuity. Hence 
X n (t) = 1 for large n. Since we have pointwise convergence on a set of (a.s.) full 
measure, we can apply the dominated convergence theorem to conclude 

lim / XJt)dt= f l B (St)dt, a.s. 
But this is the desired result, since 

r T n pT 



/ X n {t)dt = J2 X m {t)dt 
Jo i=1 Jo 

n „ 

=J2 x ni (t)dt 



n 1 n 

i=l i=l 



□ 



Theorem [4] suggests the approximation 
(28) e- rT V(F - A)+ w ne- rT E ^ - jf l B (5 t )d^ 

for the price of the structure floor (22). It is obtained from replacing F by F/n in 
the relation 

E(nF — A) + ~ nE - ^ l B (S t )cftJ , 

which follows from Theorem |lj On the right hand side of ( [28] ) we recognize the 
price of a put on the occupation time of {St), also called a corridor option. Fusai [6] 
studied such options in the Black-Scholes model. In particular, his Theorem 1 gives 
an expression for the characteristic function of ls(St)dt. Since the formula is 
rather involved, we do not reproduce it here. Section 4 of Fusai [6] explains how 
to compute the corridor option price from the characteristic function by numerical 
Laplace inversion. 

This approximation holds for period lengths tending to zero. One could also 
let the number of coupons tend to infinity for a fixed period length P, so that 
maturity increases linearly with n. The dependence of the random variables Ci 
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and Cj decreases for large \i — j\, and so we conjecture a central limit theorem, i.e., 
that 

A - E[A] 

converges in law to a standard normal random variable as n — > oo. Note that 
E[A] = £" =1 E[Cy and E[A 2 ] = E[A] + 2 J2 iKj MQCj] can be easily computed 
from Proposition [l] respectively Theorem [2] The structure floor ( 22 ) could then be 
approximately valuated by a Bachelier-type put price formula. We were not able, 
though, to verify any of the mixing conditions [4 that could lead to a central limit 
result. This is therefore left for future research. 
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